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Abstract. The canonical Schmidt decomposition of quantum states is 
discussed and its implementation to the Quantum Computation Simu- 
lator is outlined. In particular, the semiorder relation in the space of 
quantum states induced by the lexicographic semiorder of the space of 
the Schmidt coefficients is discussed. The appropriate sorting algorithms 
on the corresponding POSETs consisting from quantum states are for- 
mulated and theirs computer implementations are being tested. 
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1 Introduction 

Let £{Ha ® Hb) be a set of quantum states of a given composite system Sa+b 
on the corresponding Hilbert space Ha <8> Hb ■ As it is well known the problem 
to decide whether a given state p € £ is entangled or not is in general NP- 
HARD [5]. Although some mathematical procedures for this purpose and also 
for deriving in a quantitative way the amount of entanglement included do exist 
the real problem with them is that they are hardly to be efficiently calculable 
[2] and [3]. 

LOCC 

In the case of pure states the operational semiorder \ip) >~ \ip) on the space 
of pure states (meaning that the state can be transformed into the vector \tp) 
by exclusive use of only LOCC class of operations) has been formulated in an 
effectively calculable way by appealing to the corresponding Schmidt decomposi- 

LOCC 

tion in Nielsen [4] . The relation >- introduces only partial order on the space 
of pure states and it is why the adaptations of standard sorting algorithms of 
the corresponding POSETs are much more involved [6] and [7] comparing to the 
case of linearly ordered sets. Several adaptations of this sorting procedures are 
being adopted and tested on the Zielona Gora Quantum Computing Simulator 
and some recent results of this kind will be presented in the present contribution, 
see also [7]. 

Another topic discussed in our contribution is an attempt to generalize the 
Nielsen result to the case of general mixed quantum states. For this purpose the 



Schmidt decomposition in the corresponding Hilbert-Schmidt space has been 
used and certain semiorder relation in the space of quantum states has been 
introduced. The effort has been made to connect the introduced semiorder with 
several notions of quantitative measures of entanglement. 

A well known distillation of entanglement procedure [1] also introduces a 
partial semiorder on the space of quantum states. However this process is hardly 
to be effectively calculable and moreover it requires to have many (infinitely 
many in fact) copies of a given unknown quantum state at hands in order to 
perform the distillation process. 

Although we have no complete proof we formulate a conjecture that the 
partial order induced by the lexicographic order of the Schmidt decomposition 
coefficients is connected to the operational meaning saying that in the state 
pi is no less entanglement contained then in the state p 2 if the state pi may 
be transformed into p 2 by means of local operations supplemented by classical 
communication only. 

2 Algorithms for sorting quantum states 

2.1 Canonical Schmidt decomposition of quantum states 

For a given finite-dimensional Hilbert space H the corresponding Hilbert-Schmidt 
space is denoted as HSiTL). Let us recall that the space HS(H) consist of all 
linear operations acting on H and equipped with the following scalar product: 



A system (Ei,i = 1, . . . , dim(H) 2 ) of linearly independent matrices on H is 
called complete orthonormal system iff (Ei\Ej)ns = 1- If moreover all E 1 , are 
hermitean the system (Ei) is called complete hermitean orthonormal system. 

Proposition 1. Let p <G £(Ha ®Hb). Then there exist: a number r > (called 
the canonical Schmidt rank of p) and a complete orthonormal system (Ef) (resp. 
{Ef ) ) in HS{H A ) (resp. HS(H B ) ) and such that 



where the numbers X a > are called (the canonical) Schmidt coefficients of p. 
If all X a are different then this decomposition is unique. 

Remark 1. A different notions of Schmidt decomposition of density matrices are 
being discussed in the literature [9]. Our Schmidt characteristics like the canon- 
ical Schmidt rank and (canonical) Schmidt coefficients and the corresponding 
orthonormal systems are in unique way connected with a given p and in princi- 
ple all the properties (separability — nonseparability for example) of p should be 
obtainable form this decomposition. For example if the corresponding E£, E^ 
in formula (2) are nonnegative and therefore hermitean then p is separable. 
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(2) 



Remark 2. It is well known [10], [11] that for separable states the sum of the 
canonical Schmidts coefficients is always less or equal to 1. This leads to the 
separability criterion known as cross norm criterion. 

Remark 3. The closed subspace of HS(H) consisting of hermitean matrices 
forms a real Hilbcrt space. Therefore if the SVD theorem extends to the real 
Hilbcrt space case then the corresponding systems in formula (2) are hermitean 
by the very construction. 

Now we formulate constructive route to the canonical Schmidt decomposi- 
tion. 

Proposition 2. Let d = dim(H) and let (Ei, i = 1, . . . , d 2 ) be a system of linear 
independent matrices on H. Then there exists operation O converting the system 
(Ei) into the orthonormal system (Fi). If the system (Ei) consists of hermitean 
matrices then 0((Ei)) is also system formed from hermitean matrices. 

Proof. The well known Gram-Schmidt orthonormalisation procedure is used as 
an example of the converting operation O. I 

Let now (F i ■ ) (resp. (Ff)) be any orthonormal system in HS(Ha) (resp.iJS'(Ws)). 
Then the system (F^ ® F?) forms a complete orthonormal system in HS(Ha <8> 
TLb)- Thus taking any p £ £(TLa ®Wb) we can decompose: 



where a tj = Tr (pF^ <g> Ff) . 

Then we apply SVD operation to the matrix C — (cj.j) yielding (like in the 
vector case) all the data for supplying the decomposition (2). In particular the 
singular values of the matrix C are equal to the squares of the Schmidt numbers 
from (2). 

2.2 Linear and partial semi-order for entanglement states 

Firstly, we present a simple algorithm to realise sorting a set of quantum states 
by using von Neumann entropy notion. We will call this algorithm a linear sort- 
ing by entropy algorithm (abbreviated as LSEA). The pseudo-code of LSEA is 
presented in Algorithm (1). 

The second presented algorithm realises sorting of entangled states using 
the Schmidt decomposition. The pseudo-code is presented in the Algorithm (3). 
The input of the Algorithm (3) is now a list V of vector states on the space 
H = Ha ® Kb. The output is divided into two parts. The first part is the 
partitioning of V: V = [Vj. , . . . , V p ] where 




(3) 



V(i) £ V, SchmidtRank(yi) = n = const; and (4) 
l<n <r 2 < ... <r p < min(dimW A , dimH B ), U t V(i) = V 



Algorithm 1 Algorithm for sorting entangled quantum states using the von 
Neumann entropy 

1: function LSEA( E : {p!,p 2 ,- ■ ■ , Pn} ) : E SORT : {p lt p 2 ,- ■ ■ , Pn} 
2: for i=l to N do 

3: P?=^H A (Pi) 

4: [<Ti, VJ] - EigenSystem(pf) 

5: En(i) = E(pf) = -£ fc A* logA« fe 

6: end for 

7: s sort _ CLASSiCALSORT({En(l),En(2) En(AT)}) 

8: return E SORT 

9: end function 




Fig. 1. General idea of sorting quantum states where the lexicographic order is used. 
In result the obtained structure represents the partial order where in the sets of quan- 
tum states with the same Schmidt rank may appear a linear chains which are non- 
comparable 



i.e. the partitioning with respect to increasing Schmidt's ranks. Additionally, we 
return the complete data of merging of each Vi . 

We need to define the corresponding semiorder relation. For a given pair of 
quantum states pi and p2 we perform the canonical Schmidt decomposition (2) 
of them obtaining the corresponding Schmidt ranks ri , r-i and the correspond- 
ing Schmidt coefficients A a (l), A a (2). Having this data we can formulate the 
following algorithm which in fact introduces the semiorder relation in the space 
of quantum states. This algorithm will be called SD-QueryOracle. 



Algorithm 2 Implementation of query oracle for sorting entangled quantum 
states 

1: function SD-QueryOracle( V : {pi,P2j ) : {(pi -< Pi),{pi < 
pi), non — comparable} 

2: for i=l to 2 do {n, A r ,(i)}} = SchmidtDecomp^) 

3: if (ri > r^) return (pi -t, p^) 

4: if (ri > ri) return (p2 -< pi) 

5: sort the sets {Ai(i)i=i,2} in non- increasing order 

6: if (ri = ra) if (V J= i,..., ri (ELi < A?(2)) ) return (pi ~< p 2 ) 

7: return non-comparable 
8: end function 



Having denned the semiorder -< we can now formulate one of the possible 
(see [12] for other versions) sorting sets of quantum states algorithm that will 
be called MergeSort type algorithm. 



Algorithm 3 Algorithm for partial ordered sorting entangled quantum states 



1: function ChainMergeSort( E : {pi, .... p N } ) : V 

2: for i=l to N do {r(i), {Ai(i), . . . , \ ri (i)}} = S(i) = SchmidtDccomp(Vi) 

3: r* — max(r(l), . . . , r(n)) 

4: for a— 1 to r* do 

5: V» = [] 

6: for i=l to N do if (r(i)=a) V(a) = [V(a), p a ] 

7: end for 

8: for a=l to r* do 

9: choice randomly p £ V(a) 

10: fl° = (p,{}) 

11: fl°(l) = p' ; fl°(l) = {} 

12: (7(a) = V(a) \ {p} 

13: while (7(a) # do 

14: choice pi (E U(ot) 

15: (7(a) = (7(a) \ {p} 

16: construct a chain decomposition C(a) — {Cf, C^, . . . , C^} of -R a 

17: end while 

18: for i— 1 to q do 

19: II: do binary search on C£* using SD-QueryOracle to find smallest element (if any) 
that dominates pi 

20: 12: do binary search on using SD-QueryOracle to find largest element (if any) 

that dominates pi 

21: end for 

22: infer all results R a of II and II into R a : 

23: r(l) = fi"(l)up 

24: R" (2) = _R=(2) u R a 

25: find a chain decomposition C a of R a 

26: construct ChainMcrgc(H a , C a ) data 

27: end for 

28: return (V = [Vi, V 2 , . . . , V r ], ChainMerge(V) = [ChainMerge(V»),i=l,. . . ,r]) 
29: end function 



2.3 Example of usage of linear sorting algorithm 

Let us consider family of Bell maximally entangled states for qubits and qudits. 
In the case of qubits, these states have the following form 

\ip) = a\00) ±/3\ll) or \tp) = a\01) ± /3\10) and \a\ 2 + |/3| 2 = 1, (5) 

and for qubits there exist exactly four such maximally entangled states. The one 
of the Bell states for qudits where d — 3 can be written similarly to the qubit 
case 

=a |00) + ai|ll)+Q! 2 |22} and |a | 2 + \ ai \ 2 + \a 2 \ 2 = 1. (6) 

In general, the set of d-level Bell maximally entangled states for two qudits can 
be expressed through the following equation: 

K> = 4=7 J2 e2mJP/d \Mj + 1) mod d). (7) 



It is possible to express equation (7) in terms of qudit gates: 



\^ d pq ) = (I d ® X d f ■ (H d ® I d ) ■ (Z d ® J d f • CNOT d • |00). (8) 

where < p < d — 1 and < q < d — 1 arc indices of one of d 2 allowed Bell state. 
The symbol / represents the identity matrix for d- level qudit, and H represents 
the Hadamard gate and Z and X are generalized Pauli's operators. 

A simple function written in the Python programming language which uses 
the QCS module to generate entangled states is depicted in Fig. 2. We use this 
function to construct entangled states for earlier prepared quantum register. The 



def make_psi(r, p, q): 


def make_psi(r, p, q): 


r. Reset () 


r. Reset () 


for i in range(0,q): 


for i in range(0,q): 


r.NotN(l) 


r.NotN(l) 


r.HadN(O) 


r.RandGateRealN(O) 


for i in range(0,p): 


for i in range(0,p): 


r.PauliZ(O) 


r.PauliZ(O) 


r.CNot(0,l) 


r.CNot(0,l) 



Fig. 2. The functions written in Python preparing the entangled Bell states for given 
register. The symbol denoted by r is an object representing the quantum register and 
p and q are indices of Bell state generated by this procedure. The left column generates 
maximally entangled state but in the right column instead of the Hadamard gate we 
use a randomly generated gate to produce a non-maximally entangled states 



function presented in the left column of Fig. 2 generates states which have always 
the same amount of entanglement. Therefore the function from Fig. 2 must be 
equipped with some additional unitary gate to modify the entanglement amount. 
In the qubit cases the additional rotation gate after Hadamard gate can be used. 
In general any unitary gate that realises the rotation through any axis may be 
used to generate Bell states with uniform distribution of entanglement. Indeed, 
the right version of function makejpsi from Fig. 2 possesses this feature. 

Using function from Fig. 2 and the appropriate computational procedure to 
calculate the von Neumann entropy it is possible to prepare a simple benchmark. 
Additionally, to obtain comparable results we prepared simple test as a script in 
Python language for quantum registers built only from qubits. The test contains 
the following computation steps: first we generate n quantum registers then for 
every register the von Neumann entropy is calculated. After these steps we sort 
the obtained list using the classical method called sorting by selection. In Fig. 3 
we present the real time necessary to perform this simple test. 



Number of registers 


Time (results in sees) 


10 


0.0008762 


100 


0.0048304 


1000 


0.1407390 


2000 


0.4907356 


4000 


1.8053686 


10000 


10.643833 



Fig. 3. The time consumed by sorting tests which use randomly generated quantum 
registers with a different amounts of entanglement in the sense of von Neumann's 
entropy 

2.4 Computational complexity analysis 

The computational complexity of Algorithm (1) is given by following equation: 

N 

) (9) 

i=i 

where N is the total number of quantum registers and di represents freedom level 
of qudit used in given m quantum register. Additionally, Ti(-) represents the 
complexity of partial trace calculation, T^-) is the complexity of calculation of 
the eigenvalues and eigenvectors and T 3 (-) is the complexity of the von Neumann 
entropy calculation. Each of mentioned complexity functions work on matrices 
and if we assume that n is the size of matrix and d is freedom level of a state 
which is given by density matrices we obtain the following relations: 

T 1 (n,d)^dn 2 , T 2 (n,e0=n 3 , T 3 (n,d)=n. (10) 

The complexity of T sort (n) depends on the algorithm used to sort the obtained 
quantum registers, the value of entropy is used to compare two registers. If we 
use one of the popular sorting methods like Heapsort with complexity given by 
0(n log (n)), the complexity of Algorithm (1) will be 

T(n) = N(dn 2 + n 3 + n) + ralog(n), (11) 

where the process of computation of eigenvalues and eigenvectors is the most 
time-consuming part of the whole process of sorting quantum registers. 

The second algorithm of sorting quantum states (Algorithm (3)) contains 
oracle routine as described by Algorithm (2). The complexity of oracle for the 
worst case when the ranks are equal is given by: 

T(n) = n 3 + n log(n) + n 2 = 0(n 3 ), (12) 

The procedure of calculation the singular value decomposition dominates the 
computational complexity of oracle routine. It is important to stress that in the 
oracle routine we also sort the Schmidt coefficients, but by using the classically 
effective algorithm. However, the SVD still dominates the complexity of SD- 
QueryOracle. 



It is known [7] that Algorithm (3) calls the query at most 0(w ■ nlogn), 
where w is the maximal width of posct containing n elements but the time of 
SVD again dominates the whole process of partial sorting of quantum states. 

3 Conclusions and further work 

Basing on the canonical Schmidt decomposition of quantum states a specific 
semiorder relation has been introduced in the space of quantum states of a given 
bipartite system. In the case of pure states the introduced semiorder relation 
possesses a very clear operational meaning as described by Nielsen [4] for the 
first time. Whether the same operational meaning can be affiliated with the 
analogous semiorder relation defined in the space of all quantum states should 
be explained. 

Additionally, some version of sorting algorithm of the arising posets, the so 
called ChainMerge sorting and basing on the particular version of query oracle 
comparing the amount of entanglement in two quantum states is presented and 
tested in the case of vector states. The following extensions of the present ma- 
terial seems to be worthwhile to perform: (a) to extend the Nielsen result [4] to 
cover the case of general quantum states, (b) to formulate several different version 
of sorting posets algorithm with special emphasis putted on their computational 
complexity, (c) to formulate different version of query oracles for comparing the 
amount of entanglement included in two general states of bipartite systems. 
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